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1.  Introduction 

Quantum  computing  originated  with  Feynman’s  observation  that  quantum  systems  are 
hard  to  simulate  on  classical  computers,  but  that  they  might  be  easier  to  simulate  if  one 
had  a  computer  which  could  be  operated  quantum  mechanically  [1].  Developments  dur¬ 
ing  the  subsequent  two  decades  have  not  only  supported  this  observation  [2—8],  but  have 
also  demonstrated  that  quantum  computers  would— if  they  existed— solve  certain  classical 
problems  like  factoring  more  efficiently  than  is  possible  on  classical  computers  running  the 
best  classical  algorithms  known  [9].  This  raises  a  natural  question  [10-12]:  Might  quan¬ 
tum  computers  help  with  the  simulation  of  classical  systems?  Or  more  specifically,  given 
the  focus  of  this  workshop,  could  quantum  computers  efficiently  simulate  fluid  dynamics? 
Without  answering  the  specific  question,  in  this  paper  I  try  to  explain  why  the  answer 
to  the  general  question  may  be  ‘yes’,  and  along  the  way  explain  some  of  the  quantum 
algorithmic  tricks  which  seem  likely  to  be  useful  in  future  investigations  of  these  questions. 

I  begin  by  (very)  rapidly  introducing  the  ideas  of  quantum  computation  in  §2.  For 
a  complete  presentation,  see  [13].  In  §3  I  describe  one  approach  to  simulating  quantum 
systems — quantum  lattice  gas  automata  (QLGA)  [14]  and  then  in  §4  explain  a  connection 
with  simulation  of  classical  systems.  The  crucial  issue  is  the  relative  complexity  of  quan¬ 
tum  and  classical  algorithms;  so  §5  contains  a  detailed  analysis  for  one  specific  problem, 
including  some  new  results  on  implementing  QLGA  algorithms  on  standard  quantum 
computers.  I  conclude  with  a  brief  discussion  in  §6. 

2.  Quantum  computers 

The  possible  states  of  a  classical  computer  are  (very  long)  bit  Strings  b\ . . .  bn  6  {0, 1}  .A 
particular  computation  proceeds  via  a  sequence  of  maps  to  new  bit  strings:  Zq . . .  bn,  . . ..  A 
fundamental  result,  which  contributed  directly  to  the  conceptual  development  of  quantum 
computation,  is  that  any  classical  computation  can  be  made  to  be  reversible  [15],  ie.,  these 
maps  can  be  chosen  to  be  permutations  on  the  space  of  states.  Quantum  computation,  can 
then  be  understood  as  a  generalization'  of  classical  mmpnmionu  Thapossible  states  of  a 
quantum  computer  are  superpositions  of  bit  strings  X  a&i...6n  l^i  •  •  *  bn)  ^  (each  C 

tensor  factor  is  called  a  quantum  bit  or  qubit),  where  X  Ki...fcJ2  =  1  so  that  the  norm- 
squared  of  each  amplitude  n  G  C  is  the  probability  that  the  state  \bi  —  bn)  is  observed 
if  the  quantum  system  is  measured  in  this  ‘computational  basis’.  A  particular  quantum 
computation  proceeds  via  a  sequence  of  unitary  maps  to  new  states  Xa6i...bnl^1  mmmbn)i 

_ This  much  is  a  generalization  of  classical  reversible  computation  since  permutations 

are  unitary  maps,  and  each  classical  state  is  an  allowed  quantum  state.  The  difference 
is  that  the  final  state  is  not  directly  available;  it  can  only  be  sampled  according  to  the 
probabilities  given  by  the  norm-squared  of  the  amplitudes. 

To  evaluate  the  computational  complexity  of  an  algorithm,  either  classical  or  quan¬ 
tum,  we  must  specify  a  set  of  elementary  operations,  the  number  of  which  used  during 
the  computation  quantifies  the  complexity.  If  we  allow  arbitrary  permutations  classically, 
or  arbitrary  unitary  transformations  quantum  mechanically,  any  state  can  be  reached  from 
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a+1 


NOT 


b+a 


C-NOT 


any  other  state  in  a  single  step — these  are  a - [T]- 

clearly  not  reasonable  models  of  compu¬ 
tation.  Instead,  consider  the  ‘gate  opera¬ 
tions’  shown  in  Fig.  1:  NOT,  C-NOT  and  C- 
C-NOT  (‘C’  abbreviates  ‘CONTROLLED’). 

Each  of  these  is  a  permutation  on  the  set 
of  bit  strings,  and  Toffoli  [16]  has  shown 
that  this  set  of  three  gates  is  universal  for 
classical  computation,  in  the  sense  that 
any  (reversible)  boolean  operation  can  be 
decomposed  as  a  sequence  of  these  gates — 
which  by  Bennett’s  result  [15]  suffices  for 
universality.  Each  of  these  gates  can  be 
extended  by  linearity  to  a  unitary  map  on 
(C2)®n,  acting  non-trivially  only  on  a  sub¬ 
set  of  1,  2,  or  3  qubits,  and  thus  can  also 
be  used  as  a  quantum  gate.  (Tradition¬ 
ally  the  quantum  NOT  gate  is  denoted  by 

ux  or  X.)  Two  other  quantum  gates  which  are  particularly  useful  act  on  single  qubits 


- € 

3 - 

c+ab 


C- C-NOT 


Figure  1.  A  universal  set  of  classical  reversible 
gates.  Inputs  are  on  the  left;  outputs  are  on  the 
right.  +  denotes  addition  mod  2. 


H 


-L(]  _\)  and  JC=(J  °), 


the  ‘Hadamard  transform’  and  phase  u  rotation,  respectively.  These  matrices  have  been 
expressed  in  the  computational  basis;  thus 


tf|0) 

tf|iy~ 


V2' 


(!o>- 


|r». 


C-NOT  and  for  u>  =  eie  with  6  =  cos  1 1  form  a  universal  set  of  gates  for  quantum 
computation  [17]. 


As  we  noted  earlier,  any  array  of  (reversible)  classical  gates  can  be  simulated  by  some 
array  of  quantum  gates.  The  remarkable  fact  is  that  in  some  cases  fewer  gates  are  required 
quantum  mechanically.  The  following  example  is  the  smallest  version  of  the  Deutsch- Jozsa 
[18]  and  Simon  [19]  problems: 

Example.  Given  a  function  /  :  (0, 1}  — >  (0, 1},  we  would  like  to  evaluate  /( 0)  XOR 
/( 1).  The  function  is  accessed  by  calls  which  have  the  effect  of  taking  classical  states  ( x ,  b ) 
to  (x,  b  ©  /(x))  where  ©  denotes  addition  mod  2.  This  is  a  reversible  operation  and  thus 
also  defines  a  unitary  transformation  on  a  pair  of  qubits:  /-C-NOT.  Classically,  this  gate 
must  be  applied  at  least  twice  (once  with  x  =  0  and  once  with  x  =  1)  in  any  algorithm 
which  outputs  /( 0)  XOR  /( 1)  correctly  with  probability  greater  than  |  (assuming  a  uniform 
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distribution  on  the  possible  functions).  Quantum  mechanically  we  can  exploit  interference 
to  do  better,  applying  the  operation  only  once.  Suppose  the  system  is  initialized  in  the 
state  |0)  0  |0).  Then  apply  the  following  sequence  of  unitary  operations: 

|0)®|0>2225£l|*>®(|0)-|l» 
x=0  Z 

I)nx)\x)®(\0)-\1)) 

x=0 

XjV—0 

=  -^|/(0)  XOR  /(1)>  0(|O)  -  '|1)). 

The  first  step  Fourier  transforms  the  ‘query’  qubit  into  an  equal  superposition  of  |0)  and 
|1),  and  initializes  the  ‘response’  qubit  into  a  state  which  will  create  the  phase  (-l)/(x) 
in  the  second  step.  The  third  step  Fourier  transforms  the  query  qubit  again,  creating  the 
interference  which  ensures  that  subsequently  measuring  the  query  qubit  outputs  /( 0)  XOR 
/(l)  correctly  with  probability  1.  The  gate  array  implementing  this  quantum  algorithm  is 
shown  in  Fig.  2;  it  includes  only  one  f-C- NOT  gate. 

|o> 

|o> 

Figure  2.  The  gate  array  implementing  the  quantum  algorithm  which  solves  the  XOR  problem 

with  a  single  call  to  the /-C-NOT  gate.  The  states  are  not  normalized. 

This  example  displays  perhaps  the  simplest  quantum  improvement  over  classical  com-, 
putation.  Other  superior  quantum  solutions  to  classical  problems  include  Deutsch  and 
Jozsa’s  balanced  function  algorithm  [18],  Bernstein  and  Vazirani’s  parity  search  algorithm 
[20],  Simon’s  period  finding  algorithm  [19],  Shor’s  factoring  algorithm  [9],  Grover’s  unstruc¬ 
tured  search  algorithm  [21],  van  Dam’s  algorithms  for  weighing  matrices  and  quadratic 
residues  [22],  and  Hunziker  and  Meyer’s  highly  structured  quantum  search  algorithms  [23]. 

3.  Quantum  lattice  gas  automata 

Although  the  example  in  §2  demonstrates  the  superior  computational  power  of  quantum 
systems  for  certain  problems,  it  seems  to  have  little  to  do  with  confirming  Feynman  s 
original  idea  that  quantum  systems  could  efficiently  simulate  other  quantum  systems  [1]. 
For  the  community  attending  this  workshop,  a  natural  place  to  look  for  such  confirmation 
is  lattice  models:  The  possible  configurations  for  each  particle  on  a  one  dimensional  lattice 
L  are  labelled  by  pairs  (x,a)  e  Lx  {±1},  where  x  is  the  position  and  a  the  velocity. 
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A  classical  lattice  gas  evolution  rule  consists  of  an  advection  stage  (x,a)  (x  +  a,  a), 
followed  by  a  scattering  stage.  Each  particle  in  a  quantum  lattice  gas  automaton  (QLGA) 
[14]  exists  in  states  which  are  superpositions  of  the  classical  states:  |  ip)  =  X.  VVal^:  a)> 
where  1  =  (ip\ip)  =  The  evolution  rule  must  be  unitary;  the  most  general 

with  the  same  form  as  the  classical  rule  is: 


^2^x,a\x,a) 


advect 


scatter 

l - > 


y%2'4’x,a\x  +  a,  a) 

^  1 Px,aSaa '  |x  -|-  Q!,  G;  ), 


where  the  scattering  matrix  is 


(  cos  s 
y  i  sin  s 


i  sin  s 
coss 


Let  U  denote  the  complete  single  timestep 
evolution  operator  for  a  single  particle, 
the  composition  of  advection  and  scatter¬ 
ing.  Fig.  3  illustrates  this  quantum  evo¬ 
lution:  at  s  =  0  it  specializes  to  the  clas¬ 
sical  deterministic  lattice  gas  rule.  The 
Ax  =  At  — ►  0  limit  of  this  discrete  time 
evolution  is  the  Dirac  equation  [14];  the  Ax2 
[24]. 


Figure  3.  The  general  evolution  rule  for  a  single 
particle  in  the  one  dimensional  QLGA. 

=  At  — >  0  limit  is  the  Schrodinger  equation 


QLGA  models  can  be  generalized  to  higher  dimensions  [24] ,  and  to  include  more  parti- 
cles  [14,25,24],  potentials  [26,27]  and  various  boundary  conditions  [28].  These  are  quantum 
models  which  we  might  try  to  simulate  classically  or  quantum  mechanically.  Figures  4-6 
show  the  results  of  classical  simulations  of  plane  waves,  wave  packets,  and  scattering  off  po¬ 
tential  steps,  respectively.  These  support  the  claim  in  the  previous  paragraph  that  QLGA 
are  discrete  models  for  quantum  particles.  Li  the  next  section  we  will  see  how  they  are  also 
relevant  to  the  question  of  simulating  classical  physical  systems  quantum  rpechanicallv. 


Figure  4.  Plane  waves  in  the  general  one  dimensional  QLGA  [26]. 
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4.  Diffusion 

The  evolution  rule  for  a  single  particle  QLGA  bears  some  resemblence  to  a  random  walk. 
More  precisely,  it  is  the  unitary  version  of  a  correlated  random  walk  [29,30]  with  6  =  7r/4 
the  analogue  of  uncorrelated.  In  Fig.  5,  for  example,  we  can  see  evolution  like  that  of  a 
biased  random  walk,  with  a  spreading  Gaussian  distribution.  And,  in  fact,  this  is  as  close 
as  possible — there  is  no  quantum  random  walk  in  the  sense  of  local,  solely  x  dependent 
unitary  evolution;  the  a  dependence  must  be  included  [31].  (This  is  not  true  for  quantum 
processes  continuous  in  time  [32],  which  can  also  be  used  for  computational  purposes  like 
the  ones  we  are  considering  here  [33,34].) 

Consequently,  there  are  differences.  Diffusion  approaches  an  equilibrium  state,  in¬ 
dependently  of  the  initial  condition  (on  connected  spaces).  Unitary  evolution  of  a  single 
particle  QLGA  cannot:  the  distance  5  between  successive  states,  defined  by  cos  6  =  {Uip\ip) 
is  constant,  so  the  evolution  of  any  state  which  is  not  an  eigenstate  of  U  with  eigenvalue 
1  does  not  converge.  Each  state  implies  a  probability  distribution  on  the  lattice,  given  by 

Pt(x)  =  prob(z;t)  =  |V'x,-i(f)|2  +  |V’*,+i(*)|2> 

so  we  can  also  ask  if  this  converges.  In  fact,  this  probability  distribution  is  constant  for 
each  eigenstate  \<pi)  of  U  [26].  But  since  there  exists  T  6  Z>o  such  that  Xf  ~  1  for  all 
eigenvalues  Aj  of  U,  and  hence  ip  ~  ip,  for  any  initial  state  such  that  Pi  ^  Po,  the 
probability  distribution  cannot  converge  either  [35]. 

Aharonov,  Ambainis,  Kempe  and  Vazirani  have  shown,  however,  that  the  time  average 
of  the  probability  distribution  does  converge  [35]:  Expand  the  initial  state  | ip)  =Yli  ai\  fa) 
in  terms  of  the  eigenvectors  of  U.  Then  Ul\ip)  =  so 

prob(rr,  or,  t)  =  |y^aiA^(a:,0!|^i)| 

=  ^  ^ a,j dj( X% A j )  (x, Oi\<Pi){<pj \x, ci). 
i,3 

Then  the  time  average  of  the  probability  is 

1  T_1  1  T_1 
—  ^  prob(x,  a;  t )  =  '^2o-iaj(x,a\(pi)(<pj\x,a)—  ^(AjA j)*. 

1  t=0  i,j  t- 0 

For  Aj  A  j,  the  interior  sum  goes  to  0  as  T  — >  oo.  This  leaves  only  the  terms  in  the 
sum  for  which  A i  =  A j,  which  are  independent  of  T.  Thus  the  time  average  converges. 
In  particular,  for  the  one  dimensional  single  particle  QLGA,  it  converges  to  the  uniform 
distribution  which  is  the  equilibrium  distribution  for  diffusion  on  one  dimensional  lattices. 
That  is,  by  measuring  the  position  at  random  times  we  can  simulate  sampling  from  the 
equilibrium  distribution  of  classical  diffusion.  Although  this  is  not  true  for  all  graphs, 
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e.g.,  the  Cayley  graph  of  the  nonabelian  group  S3  [35],  we  have  analyzed  one  example  of 
discrete  quantum  simulation  of  a  classical  physical  process. 


5.  Computational  complexity 

To  be  truely  useful,  the  quantum  computation  should  be  more  efficient  than  the  corre¬ 
sponding  classical  computation.  Classically,  0(N 2)  steps  of  a  random  walk  are  required 
to  approximate  the  equilibrium  distribution  on  a  lattice  of  size  N.  Aharonov  et  al.  have 
shown  that  only  0(N  log  N)  steps  of  the  single  particle  QLGA  are  required  for  equally 
accurate  sampling  [35],  and  more  detailed  calculations  by  Ambainis,  Bach,  Nayak,  Vish- 
wanath  and  Watrous  [36]  show  that  0(N )  steps  suffice.  The  proofs  of  these  results  de¬ 
pend  on  careful  estimates  about,  for  example,  the  distribution  of  eigenvalues  of  U  and 
are  somewhat  involved.  Simple  simulations,  however,  provide  a  heuristic  explanation  for 
this  quantum  improvement.  Fig.  7  shows  the  evolution  of  a  QLGA  particle  initialized 
at  \rp)  =  ([0,-1)  -I-  |0,+l))/>/2.  Notice  that  the  peaks  of  the  probability  distribution- 
indicated  by  the  darkest  squares  in  the  plot-spread  approximately  linearly  in  time.  This 
is  the  origin  of  O(N)  number  of  steps  required  for  the  QLGA  to  sample  the  equilibrium 
distribution. 

Thus  the  QLGA  provides  a  quadratic  improvement  in  the  number  of  steps  required 
compared  to  a  classical  random  walk.  To  verify  a  computational  improvement,  however, 
the  QLGA  must  not  require  much  more  computation  per  step.  Each  step  of  the  random 
walk  requires  a  coin-flip,  i.e.,  generation  of  a  random  number.  Inasmuch  as  this  is  possible 
classically,  it  requires  constant  time,  independent  of  N .  Each  step  also  requires  the  addition 
of  the  ±1  result  of  the  coin  flip  to  the  current  position.  Since  the  latter  is  a  log  N  bit  integer, 
this  requires  O(logA)  elementary  operations.  Thus  the  total  number  of  computational 
steps  to  simulate  one  run  of  the  random  walk  is  0(N2  log  N).  We  can  compare  this  with 
the  computation  required  to  compute  the  whole  probability  distribution  by  evaluating  the 
Markov  process.  Although  this  could  be  computed  by  matrix  multiplication  at  each  step, 
the  locality  of  the  process  means  that  we  need  only  compiite 

prob(:r;  t  +  1)  =  ^prob(a:  -  1;  t)  +  iprob(a;  +  1;  t ) 

for  each  lattice  point.  Since  this  has  a  constant  computational  cost  per  lattice  point, 
evaluating  the  whole  probability  distribution  for  0(N2)  steps  takes  0(N3)  elementary 
operations. 

For  the  QLGA  we  cannot  run  single  trajectories  since  that  would  miss  any  interference 
between  trajectories.  Thus  classical  simulation  of  the  QLGA  must  be  like  the  Markov 
process  calculation  of  the  whole  probability  distribution  of  the  random  walk.  Again,  by 
locality,  each  step  requires  constant  computation  per  lattice  point.  Thus  evolution  of  the 
whole  state  for  O(N)  steps  takes  0(N2)  elementary  operations.  Taking  the  time  average  to 
approximate  the  diffusive  equilibrium  distribution  requires  another  O(fV)  factor,  hence  the 
same  0(N 3)  elementary  operations  we  found  in  the  previous  paragraph.  Unsurprisingly, 
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Figure  7.  Evolution  of  a  single  particle  QLGA  initialized  at  x  =  0  in  an  equal  superposition 
of  left  and  right  moving  states  [14].  Each  square  corresponds  to  one  of  the  classical  states;  thus 
each  lattice  site  is  represented  by  two  squares,  one  for  each  allowed  velocity.  The  boundary 
conditions  are  periodic  and  time  runs  upward. 
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therefore,  we  have  not  discovered  a  faster  classical  algorithm  for  simulating  the  equilibrium 
distribution  of  diffusion. 


Thus,  to  realize  a  quantum  mechani¬ 
cal  improvement  we  must  be  able  to  com¬ 
pute  a  single  step  U  of  the  QLGA  evo¬ 
lution  with  fewer  than  0{N )  elementary 
operations  on  a  quantum  computer. 
Schematically,  U  is  implemented  as  shown 
in  Fig.  8,  acting  on  log  N  qubits  which 
encode  the  position  on  the  lattice  and  a 
single  qubit  which  encodes  the  velocity. 


U 

H 

-a — 

Figure  8.  The  high  level  circuit  for  U.  There 
are  log  N  qubits  labelling  the  position,  and  a  last 
qubit  encoding  the  velocity. 

Advection  is  a  shift  operation,  i.e.,  mul- 


The  scattering  operation  S  acts  on  the 
last  qubit  and  thus  has  constant  cost  per  step, 
tiplication  by  a  matrix  which  is  diagonalized  by  the  discrete  Fourier  transform  FN.  The 

left  shift  |s)  • — »•  |ac  —  1)  is 


/°  1 

\ 

/I 

0  1 

U 

l  ^ 

0  1 

=  Fn 

UJ 2 

Vi  0/  V 


\ 


uN~x  / 


(1) 


Here  u>  =  e1'KilN  and  the  right  shift  diagonalizes  to  the  same  matrix  using  rather  than 
Fn  .  Assuming  N  =  2n,  the  diagonal  matrix  can  be  implemented  by  n  =  log  N  single  qubit 
phase  rotations:  R <8>  Rw2  ®  Ru4  ®  '  ‘  ’  ®  RwN /2  • 

Naively,  the  Fourier  transform  Fn  re¬ 
quires  0(Ar2)  .operations.  The  classical 
FFT  improves  this  to  0(N  log  N).  And 
perhaps  the  most  fundamental  result  in 
quantum  computation  is  that  Fn  can  be 
implemented  with  0(log2  N)  elementary 
quantum  operations  [9,37].  Thus  we  can 
implement  the  advection  step  with  an  a- 
controlled-shift  operation  |x,  a)  *— >  \x  + 

<*,  a),  where  the  shift  (in  each  direction) 
is  decomposed  as  in  (1).  Slightly  less  effi¬ 
ciently,  but  with  the  same  0(log2  N)  com¬ 
plexity,  we  can  implement  it  by  acting  by  a  left  shift,  and  then  acting  by  a  two  step  right 
shift,  conditioned  on  a  =  1;  this  circuit  is  shown  in  Fig.  9.  Thus  on  a  quantum  computer 
we  can  sample  the  diffusion  equilibrium  with  0(N  log2  N)  elementary  operations.  This  is 
an  improvement  by  0{N /  log  iV)  over  the  classical  random  walk  algorithm. 


Figure  9.  A  low  level  implementation  of  U .  The 
phase  rotations  are  as  in  the  text,  namely  R w  ® 
Ru2  ®Ru 4  in  each  of  the  three  timesteps  at  which 
they  are  applied. 


10 


Quantum  computing  classical  physics  David  A.  Meyer 

6.  Discussion 

We  have  seen  that  a  classical  physics  problem — sampling  from  the  equilibrium  distribution 
of  a  diffusion  process — is  solvable  more  efficiently  with  a  quantum  computer  than  with 
a  classical  one,  in  the  sense  that  the  QLGA  algorithm  outperforms  the  random  walk 
algorithm.  The  solution  utilizes  two  of  the  fundamental  quantum  speedups:  First,  the 
quadratic  improvement  in  the  number  of  steps  necessary  for  a  QLGA  simulation  compared 
to  a  random  walk  simulation  is  reminiscent  of  the  quadratic  improvement  of  Grover’s 
quantum  search  algorithm — which  he,  in  fact,  describes  as  diffusion  [21].  Second,  the 
exponential  improvement  of  the  quantum  Fourier  transform  over  the  FFT  [9,37]  provides 
a  speedup  in  the  advection  step.  Perhaps  the  simple  problem  considered  here  will  inspire 
application  of  these  techniques— or  new  ones — to  speed  up  computation  of  other  classical 
systems  with  quantum  algorithms.  For  example,  a  natural  generalization  would  be  to 
analyze  diffusion  in  a  potential,  for  which  a  QLGA  simulation  has  already  been  shown 
to  evolve  the  mean  of  the  distribution  quadratically  faster  than  does  a  classical  (biased) 
random  walk,  in  certain  cases  [30].  Successful  development  of  such  quantum  algorithms 
should  provide  additional  impetus  to  efforts  towards  building  a  quantum  computer  [38]. 
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